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We use the Gorkov formulation of the Dirac-Hartree-Fock-Bogoliubov approximation to nuclear 
pairing to study the x So nucleon-nucleon correlations in nuclear matter. We find the short-range 
correlations of the 1 So pairing fields to be almost identical to those of the two-nucleon virtual state. 

We obtain mutually consistent results for the pairing fields, using several different sets of effective 
interaction parameters, when we demand that each of these sets places the virtual-state pole at its 
physical location. 


The non-relativistic BCS and Hartree-Fock-Bogoliubov (HFB) approximations have been used with success to 
study pairing in nuclear physics, both years ago 01 and still today [|] [h|. Pairing approximations provide a 
simple means of extending the independent-particle approximation to one which includes the effects of the binding 
energy and short-range correlations associated with bound pairs of nucleons in the nuclear medium. The long-range 
correlations associated with collective phenomena in nuclei are similarly described using the non-relativistic random- 
phase approximation (RPA). 

On the other hand, it is well known that non-relativistic independent-particle approaches using realistic two-body 
interactions have difficulties in accounting for basic phenomena, such as the spin-orbit part of the nucleon-nucleus 
interaction and the saturation properties of nuclear matter. These properties can be fairly easily described in a 
relativistic formulation, in which effective mesons are exchanged between Dirac nucleons. The success of the relativistic 
mean field approach, initially developed by Walecka and collaborators and later by many others §| , invited its 
extension to approximations that could take into account the residual correlations between nucleons. 

Consistent relativistic formulations of both the RPA and the HFB 25-28) have thus been developed. Of the 
latter, the authors of Ref. (^5| derived self-consistent equations for the components of the relativistic pairing field but 
performed no calculations. The results obtained in Refs. |2(| and ^7j, using a zero-range model and a non-relativistic 
reduction of the pairing equations, respectively, are generally much larger than those obtained in non-relativistic 
calculations, although both of the relativistic calculations used interaction parameters that were fit to yield the 
saturation point of nuclear matter. Such results would indicate that other aspects of the interaction are involved in 
pairing than those necessary to describe the bulk characteristics of nuclei and of nuclear matter. The crudeness of 
the calculations, however, do not permit a firm conclusion to be reached. The results obtained in the fully relativistic 
finite-range calculations of Ref. |28| are very similar and do permit such a conclusion. Here we will show that, 
in addition to describing the bulk properties of nuclear matter, a nucleon-nucleon interaction must provide a good 
description of the two-nucleon x So virtual state, if it is to furnish a good description of pairing in nuclear matter. 

The HFB approximation introduces short-range two-nucleon correlations through the pairing fields and their mo¬ 
mentum dependence. These correlations have been investigated in recent years in non-relativistic models of nuclear 
pairing. jTTj . [T~i| [T7i| In particular, the close association between the short-range correlations of the two-nucleon 1 So 
virtual state and those of the 1 So pairing fields was investigated in Ref. 0- We investigate the same association in 
the relativistic model and reach a similar but stronger conclusion: the short-range correlations of the pairing field are 
almost identical to those of the two-nucleon virtual state. For the relativistic interactions we have used, the short- 
range repulsion and medium-range attraction resulting from the exchange of effective mesons provide a description of 
the short-range correlations in very good agreement with those of the two-nucleon virtual state, deduced from realistic 
nucleon-nucleon potential models. By requiring that the interactions also reproduce the position of the virtual state 
pole or, equivalently, the two-nucleon singlet scattering length, we will show that the effective relativistic interactions 
can provide a consistent description of pairing, both among themselves and with the non-relativistic calculations. 


PAIRING AND TWO-NUCLEON CORRELATIONS 

In a two-particle system, bound-state correlations can be roughly classified as asymptotic ones and short-range 
ones. The asymptotic ones are determined principally by the binding energy, while the short-range ones depend on 
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the high-momentum components of the wave function. This can be seen by examining the manner in which a bound 
pair appears in the two-body T-matrix, T(E). The T-matrix satisfies the integral equation 

T{E)=V + VG 0 (E)T(E), (1) 

where V is the two-body interaction and Gq(E) the free two-body propagator. A bound state appears in the T-matrix 
as a pole at a negative value of the energy, E = — e b , where e b is the binding energy of the pair. We then have 

T(E) = d-=^—dP+T c {E), (2) 

with T C (E) the continuum (positive energy) part of the T-matrix and d the bound-state vertex function. Substitution 
of the latter expression into the integral equation immediately yields an equation for the vertex function, 

d = VG 0 (-e b )d. (3) 

The vertex function describes the momentum dependence (spatial dependence) of the bound state and is closely 
related to its wavefunction. By rewriting the equation for the vertex function as 

d = (-£6 - H 0 ) G 0 (—Q>) d = V G 0 (-e b ) d , (4) 

we can identify the bound-state wave function as ip = Go(—e b )d. In this last expression, we can see the rough 
division of the correlations into asymptotic ones determined by the singularity of the Green’s function, Gq(— e b ), and 
short-range ones contained in the high momentum components of the vertex function, d. 

The two-nucleon system in the vacuum has a bound state in the isospin one, 3 Si- 3 Di channel — the deuteron. 
In the isospin zero, 3 So channel, where we will study pairing, the two-nucleon system in the vacuum has no bound 
state. It does have a virtual state, however, with k v ~ —0.05* fm _1 and e v = k%/M rts 140 keV, corresponding 
to a two-nucleon singlet scattering length do = l/|fc.„| « 23 frn. An expansion similar to that of Eq. (|j) can still 
be performed to extract the contribution to the T-matrix of the virtual state. Because the energy of this state is 
extremely small, its contribution to the T-matrix dominates the low-energy scattering and essentially determines the 
short-range correlations in the 3 So channel. The short-range correlations determined by the momentum dependence 
of the vertex function of the virtual state, d v (q ), are thus to good approximation given by the momentum dependence 
of the half-on-shell 3 So T-matrix at zero energy, fi O (g,0;0). We have from the expression analogous to Eq. d). 

tio(q, 0; 0) =< q\T w (E = 0)|A; = 0 >«< q\d v > — < d„|0 > ■ (5) 
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Note that the asymptotic properties of the virtual state also play a role here, by determining the magnitude of the 
contribution of the virtual state to the T-matrix through the factor l/e„. The short-range correlations contained in 
the g-dependence can be extracted unambiguously, however, by normalizing both sides of the equation to their values 
at q = 0. 

To look for general two-particle correlations in nuclear matter, one could study the Brueckner C/-matrix. In this 
generalization of the two-particle T-matrix of Eq. (Q), the propagator Gq{E) is now a many-body operator that, 
beside describing the two-body propagation, must take into account the effects of Pauli blocking and of interaction 
with the nuclear medium. To look for bound-state correlations, one would examine the poles and respective vertex 
functions of the t/-matrix as advocated years ago by Emery . Such an approach has been used by various authors 
p9| j3l[ . We will use instead the HFB approximation. It too has been used to study bound-state correlations in 
nuclear matter |ll|,|l^,|l6| and to even estimate deuteron production in heavy-ion collisions 0' 

We can lend force to our use of the HFB approximation to study bound-state correlations by analyzing qualitatively 
the manner in which the one-body HFB approach takes into account the two-body pairing bound state. We begin by 
observing that a bound-state vertex function d of the T- or Q -matrix can be considered as an operator that converts a 
two-body bound state into two single particles. We could instead consider this a one-body operator by (1) neglecting 
the particle number of the bound state and (2) considering one of the outgoing particles as an entering hole. The 
effect of the vertex function would then be the conversion of a hole to a particle. This is just what the pairing field, 
A, of the HFB approximation does. Analogously, we can associate the adjoint vertex function, d', which converts two 
single particles to the two-body bound state, with the conjugate pairing field, A, which converts a particle to a hole. 

To obtain a complete one-body description, we must take into account the propagation of particles and holes as 
well as the conversion of one to the other. A HFB formalism that succeeds in unifying these ingredients simply and 
clearly is the Gorkov one p^]. Although infrequently used in non-relativistic studies of pairing, the formalism has 
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served as the basis for various relativistic studies mom- This is due, at least in part, to its natural expression in 
terms of the propagator language common to field-theoretical approaches. 

The Gorkov formulation of pairing extends the usual particle propagator G{x — x') to one of the form 

/ G(x — x') F{x — x')\ 

\ F(x — x') G(x — x') J ’ 

in which G(x — x') describes the propagation of holes in the medium and the anomalous propagator F(x — x') and 
its conjugate F(x — x') describe the conversion of holes to particles and particles to holes, respectively. Inverting the 
reasoning that lead us from the vertex function d to the pairing field A, we can interpret the anomalous propagator 
and its conjugate as terms describing the overlap between the two-particle bound state and the two single particles. 
We thus expect the anomalous propagators to contain information about the relative motion of the two particles in 
the bound pair. 


THE GORKOV FORMALISM 


We sketch here the development of a Dirac version of Gorkov’s self-consistent particle-hole propagator 
this, we begin with the following ansatz to the effective single-particle Lagrangian, 


28 . To do 


dtL e ff = / drxdrx' {il>(x)(i r y l/ d 1 ' — M)ip(x)8{x — x') + iup(x)~ioip(x)8(x — x') 


—ip(x)T,(x — x')ip{x') + ^if>(x) A(x — x')iP t{x') + i ip T {x)A{x — x')ip{x') 


( 6 ) 


where E is the usual self-energy, A and A are the pairing fields and fi is a chemical potential, which will be used to 
constrain the average baryon density. 

The hole wave function, ipT, is defined as 

ipT = Bip ip T = ip T , (7) 

where ip T denotes the transpose of the wave function ip, and the matrix B = T 2 ® 75 C, in which the Pauli matrix 72 
acts in the isospin space and C is the charge conjugation matrix. We note that, up to a factor of 70 and a phase, the 
isospin doublet ipT is the time-reverse of ip. 

The requirement that the effective Lagrangian be Hermitian yields the following conditions on the self-energy and 
pairing fields, 


E(cc) = 7 0 S t (-a;) 7 o and A{x) = 7 0 A t (-a:) 7 o . 

The requirement of invariance under transposition of the pairing terms yields the additional conditions 
A(x) = — B t A t (—x)B^ and A (a;) = —BA T {—x)B*. 


( 8 ) 

(9) 


These constraints are important in limiting the possible structure of the self-energy and pairing fields. 

Making use of the relation between ipT and ip, we can manipulate the effective one-particle Lagrangian into a matrix 
form, 


/ 


dtL 



d^xd^x' (^ijj(x),'ip T (x) > j x 


( 10 ) 


where 


( 




M + ii^o)5{x — x ') — S(x — x') 

A(a; — x') {i’lvd 1 ' + M 


A(x-x') 

^7o)<5(z - x') + T, t {x - x')) J V ) ’ 


Y>t{x) = BE T (-a;)B t . (11) 

The extended vector wave function is that of the quasi-particle. Its equation of motion is evident and could be written 
down immediately. 
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Due to the translational invariance of nuclear matter, the equation in momentum-space is diagonal in the wave 
number. For the momentum-space representation of the corresponding generalized Feynman propagator, the equation 
of motion takes the form 


( 'ivk v - M_— £(fc) + /rqo A(fc) 

A [k) j„k 1 ' + M + S T (fc) - /rqo 


S F {k) = 1. 


( 12 ) 


We note that the momentum-dependent unitary transformation that diagonalizes the matrix operator is the Dirac 
version of the Bogoliubov-Valantin transformation. 

To make contact between the effective quasi-particle Lagrangian and an interacting one, we assume that the 
meson+interaction terms in the latter have been reduced to four-fermion terms of the following form, 

J dtLi = — J d 4 xd l x'ip(x)r a (x)'ip(x)D a P(x — x')i!)(x')Yp{x')%l)(x '), (13) 


where Tq^x) and F^x') are vertex functions, D a @(x — x') is the meson propagator, and a and /3 represent any indices 
necessary for the correct description of the meson propagation and coupling. 

We can obtain the mean field contribution of this interaction term by replacing each of the possible pairs of fermion 
fields by its ground-state expectation value, 

J dt{Li) e ff = i J d 4 xd 4 x'D af3 (x — x'){ 2ip(x)T a (x)ip(x) < ^(x , )r / 3 (x')'*/>(x') > 

+ 2if>(x)r a (x) < ip(x)ip(x') > Yp(x')4>(x') 

- V>(a:)r a (a:) < i/j{x)i/j t (x') > r^(x , )4> T (x') 

- i/j T {x)rl(x) < -lp T (x)ip(x') > r 0 (x')ip{x')} , 

(14) 


where, by < ... >, we mean the time-ordered ground-state expectation value, < 0 | T(...) | 0 >. 

We note that the first term in this expression is a Hartree one and the second a Fock exchange one, while the last 
two, after using the definition of ipT to replace the transposed ip's, can be recognized as pairing terms. Comparing 
these mean field contributions to those of the effective quasi-particle Lagrangian, we can express the self-energy and 
pairing fields in terms of the two-fermion ground-state expectation values as 

£(a: — x') = —S(x — x , )F Q ,(a;) J d 4 x"D a ^(x — x") < ip{x")Tp{x")tp{x") > (15) 

— T a (x)D al3 (x — x') < ip(x)ip(x') > Tp(x '), 


and 


A(x — x') = —T a (x)D a/3 (x — x') < ip(x)il> T (x') > , (16) 

where the equation for A can be obtained using the Hermiticity condition of Eq. (||). These expressions become 
self-consistency equations when we evaluate the expectation values by using their relationship to the quasi-particle 
propagator, 


iSp{x — x') 


G(x — x') F{x — x') 
F(x — x') G{x — x') 


ip(x) 

ip T (x) 


(ip{x') ,4 > t {x')) 



(17) 


which is itself a function of the mean fields. 

The equations in momentum space are obtained by Fourier transforming the above expressions, giving 

£(fc) = F a (0)D^(0) J ^jTr [r^(0)G(g)] e^° 0+ (18) 

-J -0^r a ( q )D^( q )G(k-q)Tp(- q ), 

and 


4 




(19) 


A{k) = ~j (04 r*(q)D a0 (q)F(k - q)BV T 0 {-q)B t. 

To complete the set of equations, we include that constraining the baryon density, 

Pb =< 7 o^ >= J ( 01 Tr I'YoGiq)] e lq ° 0+ . (20) 

Solving the self-consistency equations in conjunction with this constraint, we obtain the nonperturbative self-energy 
and pairing fields. 

Here, we restrict our attention to pairing in symmetric nuclear matter. The Hermiticity and transposition invari¬ 
ance conditions of Eqs. (^) and (|^), as well as the requirements of invariance under Lorentz and parity transformations 
reduce the possible form of the self-energy held to 

£(fc) = E s (fc) - 70 S 0 (fc) +7 • fcSy(fc), (21) 

while we take the form of the pairing held to be 

A (k) = A (k) = (A s {k) - 70 A 0 (fc) - 1707 • kA T (k))f ■ n. (22) 

Both the self-energy and pairing terms have Lorentz scalar and time-like vector components. To be consistent with 
our assumption of symmetric nuclear matter, the self-energy must be an isoscalar. The transposition invariance 
condition of Eq. (|tj|) forces the scalar pairing held to be an isovector. We have simplihed the form of the pairing held 
by assuming that it can be taken to be real and that the isospin dependence can be isolated in an overall factor of 
t • h, where the unit vector n is arbitary. The special case t ■ h = t% corresponds to the standard one of proton-proton 
and neutron-neutron pairing. Although we have not studied more elaborate isospin dependences, we have examined 
complex solutions to the pairing equations and found them to differ by only an overall phase from the solutions 
restrained to be real. 

After substituting the simplihed expressions for the mean helds into the propagator and the self-consistency equa¬ 
tions, the latter can be reduced to coupled equations for the components of the mean helds and the chemical potential, 
p. The self-consistency equations contain contributions from both the negative-energy and the positive-energy states 
(from the Dirac sea and the Fermi sea). To avoid the complications of renormalization, we have discarded the poles 
of the HFB propagator corresponding to negative-energy states. This procedure is not equivalent to the neglect of 
the HF negative-energy states performed in Ref. ^7|, nor are the differences small, as one might hrst expect. The 
contribution to the HFB states of the negative-energy HF states reduces the magnitude of the pairing helds, much 
as the contribution to the HF states of the negative-energy free states reduces the attractive scalar component of the 
HF mean held. More details can be found in Ref. |^8), where we compare the two approximations and show that 
discarding the HFB negative-energy states yields much more reasonable results. 


NUMERICAL RESULTS AND DISCUSSION 

We have performed calculations of 1 So pairing in symmetric nuclear matter for various sets of interaction param¬ 
eters (meson-nucleon coupling constants and meson masses 0®)- We introduced, as an additional parameter, 
a momentum cutoff at \k\ = A (in the nuclear matter rest frame), which limits the momentum integrations in the 
self-consistency equations. Such a cutoff could be considered a crude approximation to the nucleon-meson vertex form 
factors that our calculations do not contain. We performed calculations for various values of the momentum cutoff. 

As the self-consistency equations for the pairing helds are proportional to the pairing helds, it is always possible 
to hnd a solution in which these helds are null. This zero pairing-held solution is just the normal HF one. In our 
calculations, we have found this to be the only self-consistent solution at sufficiently high baryon density. At densities 
lower than about two thirds of the saturation density, ps < 2pso/3, we also hnd a non-trivial HFB solution. (We 
will not consider the region of exponentially small pairing helds discussed in Ref. fl5||.) 

We display, in Fig. |lj, the two principal components of the pairing held, As(&f) and A 0 (kp), evaluated at the 
Fermi momentum, obtained for two different values of the momentum cutoff A using the a-co model parameters of 
Ref. fyj. The third component of the pairing held, At, is three orders of magnitude smaller than As and Ao and 
is not shown. The components of the pairing held are plotted as functions of the Fermi momentum, where the latter 
is defined through its standard relation to the baryon density, that is, pp = 'ykp/Gir 2 where 7 is 2 for nuclear matter 
and 1 for neutron matter. We observe that the magnitudes of the components of the pairing held, represented here 
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by their values at the Fermi momentum, Ao iS (/cf)) depend strongly on the Fermi momentum (and thus the baryon 
density). In both of the cases shown, which are typical ones, the fields increase rapidly at small densities, reach a 
maximum at about 1/4 the saturation density, and fall back to zero before the saturation density is reached. We 
observe that the magnitude of the components of the pairing field also clearly depend on the value of the momentum 
cutoff. The form of the functional dependence of the fields on the Fermi momentum, however, is fairly insensitive to 
the value of the cutoff. 

Diagonalizing the equation of motion, Eq. (|l2|), in the particle-particle and hole-hole subspaces — equivalent to 
expanding in the HF basis in these subspaces — reduces the pairing fields to effective ones for the positive- and 
negative-energy states and to a coupling term which is small for low momenta (but increases with the momentum) 
[|28| . The effective pairing field for the positive-energy states, which we call the gap function, is 


A g (fc) = 


M*(k) 

Bt 


A 0 (fc) - A s(k) + i 


k * | k 
Ei 


Ar(fc), 


(23) 


where k* = (1 + Ey(A;))|fc|, M * = M + E 0 (fc), and E* = Vfc * 2 + M* 2 . This is the quantity whose role is closest 
to that of the non-relativistic pairing field. Like the non-relativistic self-energy, it too is the difference between two 
larger relativistic quantities. This is evident in Fig. |l], in which the gap functions for the two different values of the 
momentum cutoff are displayed together with the components of the pairing field. We observe that different values 
of the momentum cut-off result in gap functions A g of different magnitude but with approximately the same form of 
functional dependence on the Fermi momentum. The gap functions are quite similar, in form, to those obtained in 
various non-relativistic calculations using realistic potentials [^|,^5 16 3(J . 

In Fig . we compare the components of the pairing field and gap functions obtained using the parameters of 
Ref. 1181 and Ref. ]lij, at a fixed value of the cutoff momentum A, again in a cr-u> model. Here, we observe that 
different values of the interaction parameters also result in pairing fields Ao and A s and gap functions A g of different 
magnitudes but, again, with roughly the same form of functional dependence on the Fermi momentum. 

As we expect the pairing fields and the gap function, like the two-nucleon vertex function, to contain only short- 
range two-nucleon correlations, we expect their dependence on the momentum to be approximately the same over a 
wide range of nuclear densities. We would expect significant changes in their momentum dependence only at densities 
high enough for the probability to become appreciable of a third nucleon to be found simultaneously with the pair 
within the effective radius of about 0.5 fm that is characteristic of the short-range correlations. In Fig. ||, we show the 
normalized momentum dependence of the pairing fields and gap function (at a fixed value of the Fermi momentum), 
A(fc)/A(0), obtained using the interaction parameters of Ref. [|19| and a momentum cutoff of A = 2 GeV/c, for several 
values of the baryon density. The curve corresponding to zero density gives the momentum dependence of the virtual 
state solution to the vacuum pairing equations, which is discussed below. Confirming our expectations, at low density, 
we find that the momentum dependence of the pairing fields and gap function is indeed almost independent of the 
density. As the density increases, differences begin to appear in the gap function but remain small almost up to the 
density at which pairing disappears. The momentum dependence of the pairing fields remains almost identical to that 
of the virtual state over the entire range of densities for which pairing occurs. 

We note that in Ref. ||l5|| a similar comparison was made of the momentum dependence of the non-relativistic gap 
function at various values of the baryon density, but for the case of neutron matter. They observed slightly larger 
variations of the normalized gap function than the ones we have found but, as their calculations were non-relativistic, 
did not observe the invariance of the pairing fields. Our results for neutron matter are almost identical to those in Fig. 
^ and are not shown here. We want to emphasize, however, that the momentum dependence of the neutron matter 
pairing fields is just as independent of the matter density as is the momentum dependence of the pairing fields in the 
case of nuclear matter. 

We would also like the short-range correlations contained in the pairing fields and gap function to be relatively 
independent of the interaction parameters, at least if we believe that these parameters provide a reasonable description 
of the actual physical situation. We can see in Fig. [| that, to a certain extent, this is the case. There we show 
the normalized momentum dependence of the pairing fields and gap function, obtained for two different values of 
the momentum cutoff using the a-ui model parameters of Ref. [ fl8| and of Ref. 0 - We see that the normalized 
momentum-dependence of the pairing fields and gap function, although similar, is clearly different for the two sets of 
interaction parameters. However, the momentum-dependence of the fields depends only weakly on the values of the 
cutoff, A, as is also the case for the other sets of interaction parameters that we have used. 

In the limit of zero density, the integral equation for the pairing field, Eq. (0), reduces to the ladder approximation 
to the Bethe-Salpeter equation for a two-particle bound state, in which the the chemical potential is identified with the 
two-body bound-state energy and antiparticle propagation in the two-nucleon center-of-mass frame has been discarded. 
For a physically reasonable interaction, this has no nontrivial solution in the channel, as the two-nucleon system 
has no bound state there. The virtual state appears as a solution when the equation is analytically extended to 
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the second energy sheet. As we have already stated, the energy of the virtual 1 So-state is extremely small, so that 
its vertex function dominates the short-range correlations in two-nucleon scattering at low energies. As the virtual 
state and the pair state at low nuclear matter densities are solutions of essentially the same equation (and contain 
essentially the same physics), we expect the short-range correlations of the pair state to be similarly dominated by 
the short-range correlations of the virtual state. We thus expect the the components of the pairing field to have 
a momentum dependence similar to that of the components of the vertex function of the virtual state. We have 
seen in Fig. |] that, at low densities, this is indeed the case. As the gap function A g has a momentum dependence 
similar to that of the non-relativistic vertex function of the virtual state d v (q) we expect, as discussed above, that it 
should also have a dependence on the momentum similar to that of the half-on-shell matrix element of the zero-energy 
non-relativistic 1 So T-matrix, tio(q, 0;0). That is, we expect that 

A g(q) _ d v (q) _ tip(q, 0; 0) 

A g (0) ~ dv(0) ~ tio(0,0; 0) 1 ’ 

In Figs. |] and [|, we have plotted, together with the normalized pairing fields, the normalized non-relativistic zero- 
energy half-on-shell 1 So T-matrix, for the case of the Bonn-B potential |h|, which includes the exchange of a, w, p, 
7 r, 77 and S mesons between nucleons. We see in Fig. |], for the interaction parameters of Ref. [Je^, that the agreement 
between the ratios of gap functions and half-on-shell T matrices is quite good. We observe in Fig. |] that the ratios 
of gap functions for the interaction parameters of Ref. [0 do not agree quite as well as those of Ref. [|id| with the 
ratios of half-on-shell T matrices, although they are still quite similar, as are the gap functions of the other sets of 
interaction parameters. We thus conclude that the momentum dependence of the effective interactions that we used 
does provide a reasonable description of the actual physical situation. 

In contrast with the almost unchanging momentum dependence of the pairing fields, we find that their magnitudes 
depends strongly on the interaction parameters and the value of the cutoff. This reflects the fact that nuclear 
pairing is the result of a delicate balance between short-range repulsion and long-range attraction. [|l3| Changing the 
interaction parameters or the momentum cutoff changes the point of this balance, modifying the magnitudes of the 
fields. However, such modifications do not appear in the magnitudes of the pairing fields alone. As we have noted, 
in the zero density limit, the equation for the pairing field, Eq. reduces to the ladder approximation to the 

Bethe-Salpeter equation for a two-particle bound state in the -^So channel. For some sets of interaction parameters 
and momentum cutoff, a bound state can indeed be found. In other cases, where no bound state exists, analytic 
continuation of the pairing equations to the second energy sheet reveals a virtual state in the channel. 

We find the magnitudes of the gap function and the pairing fields at low densities to be strongly correlated with 
the location in the complex momentum plane of the virtual state or bound state in the vacuum that is produced by 
each set of interaction parameters + cutoff. We display these correlations in Figs. [ 5 ] and |(j|, where comparison is made 
between several different parameter sets with varying numbers of mesons. Calculations were performed including the 
d, u>, p and 7T mesons and also including just the cr-w mesons, as the latter dominate the gross features of the pairing 
fields and gap function. Calculations were also performed using a relativistic zero-range interaction, in which the 
a and u> mesons are taken to be infinitely heavy and the coupling constants are defined so as to yield the correct 
nuclear-matter saturation point. For each parameter set, the momentum cutoff A was varied with all other parameters 
kept fixed. For each value of the momentum cutoff, we determined the pairing gap at £^=0.5 fm -1 , about 1/8 of the 
saturation density, and obtained the bound or virtual state energy in the vacuum. 

We observe in Figs. |5| and |(i| that the value of the free two-nucleon bound or virtual state energy essentially 
determines the magnitudes of the gap function A g (kp) and the pairing fields Ao iS (/cf) at kp=0.5 fm -1 . The a-u> 
results agree well among themselves, with their values for each of the pairing fields and gap function all lying within 
a narrow band. Although the inclusion of other mesons increases the size of the bands, their widths still remain 
small relative to the magnitudes of the fields. In Fig. ^j, we see that even the zero-range model yields results for the 
gap function in good agreement with those of the other calculations. However, it cannot reproduce the values of the 
components of the pairing field, Ao iS (/cf), obtained with the finite-range models, as is apparent in Fig. |j. 

The correlations between the components of the pairing field, Ao, s (£\f), and the position of the virtual state, shown 
in Fig. yield two distinct lines. In each case, one of these lines contains most of the finite-range results while the 
other contains the zero-range calculations. Several of the finite-range points lie between the two lines. The finite- 
range calculations that fall close to the zero-range line all correspond to extremely low values of the momentum cutoff, 
A < 500 MeV/c. They are thus similar to the momentum-independent zero-range results in the sense that the pairing 
fields (but not the gap function) in these cases vary little as a function of the momentum. (Better said, they have 
very little room in which to vary.) This is clear from the leftmost portions of Figs. and [|, in which the momentum 
dependence expected of the pairing fields for values of the cutoff A < 500 MeV/c is displayed in the range of values 
k < 2.5 fm -1 . 

We can understand better the correlation between the magnitudes of the pairing fields and the location of the 
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virtual state if we examine the coherence length of the pairing fields. The coherence length £ of a bound pair is 
defined as its root mean square radius. In terms of the momentum-space wave function of the pair, y(fc), we have 

f 2 = / // A*'®*© • (25) 

In the Bethe-Salpeter equation, the pair wave function can be expressed in terms of the residue of the two-particle 
propagator at its pole. By analogy, we argue that here we can obtain the (un-normalized) pair wave function as the 
residue of the anomalous propagator F(k) at its pole, 

X(k) = ^ f ckoF(k,u)e iu0+ , (26) 

where we retain only the contribution of the positive-energy pole, to be consistent with the vacuum truncation 
discussed earlier. In Fig. [t[ we plot the inverse of the coherence length, l/£, at fcF=0.5 fm -1 , as a function of the 
position of virtual or bound state in the vacuum, just as was done in Figs. |E| and Here too, we observe a similar 
strong correlation. What is most important at the moment, however, are the typical values obtained for the coherence 
length. We observe that, in the range of calculations shown, these never fall below 3 fm. Near the physical position of 
the virtual state, the value of the coherence length lies between about 4 to 10 fm. As will be seen shortly, the Fermi 
momentum at which we have shown the various correlations is close to that for which the coherence length of the pair 
attains its minimum value. Thus we can claim that the coherence length of the physical bound pair never becomes 
smaller than about 4 fm. 

We can thus justify the correlations observed with the following argument: For the parameter sets of physical 
interest, the size of a bound pair at low density is much greater than the range of the interaction. As the short-range 
pair correlations are about the same for all sets of finite-range interaction parameters, it is then a single parameter, the 
net attraction, that determines the differences in the binding energy of the pair obtained with each of the parameter 
sets. It suffices to fix the size (or energy) of the pair at one value of the density to obtain its trend over the entire 
range of values for which the details of the interaction can be neglected. We choose to fix the net attraction between 
pairs by fixing the position in the complex momentum plane of the virtual state in the vacuum. As we have seen, 
this is equivalent to fixing the two-nucleon singlet scattering length. 

We thus fix a momentum cutoff for each parameter set so that it places the 1 So virtual state at its physical location 
of k v « — 0.05i fm -1 . This is impossible for the a-ui interaction parameters of Ref. |n|, for which ik v > 0.08 fm -1 for 
A > 500 MeV/c. In this case, we fix the virtual state at k v ~ — 0.08i fm -1 , its closest point to the physical location. 
The pairing fields that result are then fairly consistent among themselves and with the non-relativistic result of Ref. 
[fio| for low values of the Fermi momentum, as can be seen in Fig. |j. Our results are also very similar to those obtained 
in many other non-relativistic calculations of both nuclear matter and neutron matter (for which the gap function is 
from 10% to 20%larger), which are not shown here f)[. |T5| . [i7i|; . Although the calculations shown are restricted to those 
using the a-oj parameter sets, the gap functions obtained using the a-u-n-p interaction parameters are quite similar. 
Even the zero-range model can describe the behavior of the Fermi-momentum gap-function at low densities, when it 
used with an appropriate cutoff A. 

We note that the good agreement between different sets of interaction parameters does not extend up to densities 
at which the magnitude of the pairing fields begins to decrease. There, differences in the details of the interplay 
between attraction and short-range repulsion become important, introducing a more elaborate dependence on the 
parameter sets. To the extent to which these details are reflected in the short-range correlations, that is to say, in the 
momentum dependence of the pairing fields, we can compare this momentum dependence with that of the half-on shell 
singlet T-matrix to choose the most “physical” parameter sets. As we have not performed relativistic calculations of 
the T-matrix, we have compared the non-relativistic calculations using the Bonn potential to our gap functions to 
conclude that the a — to and a — u — it — p parameter sets of Ref. [|l9| are equally good and are the most “physical” 
of those we have studied. The half-on-shell singlet T-matrix obtained using the Bonn potential and the gap function 
obtained using the <j — u parameters of Ref. [|h| are shown in Figs. |] and |j. 

We show in Fig. |]the coherence lengths £ obtained using the physical values of the momentum cutoff that resulted 
in the gap functions of Fig. [8| We observe that the coherence length drops rapidly from infinity, at low densities, 
to reach a minimum of about 5 fm at about 1/10 the density of saturated nuclear matter. It rises again rapidly, 
becoming infinite at the density at which pairing disappears. The minimum coherence length is fairly independent 
of the interaction parameters and is about 5 or 6 fm in all cases studied. The coherence length in neutron matter 
displays a very similar dependence on the Fermi momentum kp but is about 10% smaller than in nuclear matter of 
the same Fermi momentum. 

We conclude by interpreting the spatial form of a bound 1 So pair in nuclear matter, when it exists, in terms of the 
magnitude and form of the pairing field, A (k). The short-range correlations of the pair, characterized by the form 



of the pairing field, A(fc)/A(0), are almost independent of model parameters, momentum cutoff and nuclear density, 
and are almost identical to those of the two-nucleon ^So virtual state in the vacuum. The asymptotic nature of the 
pair, is roughly determined by the magnitude of the pairing fields at the Fermi momentum, which are strongly varying 
functions of the nuclear matter density. Based on the variations observed in the pair coherence length, we estimate 
that a 1 11 S 0 pair in nuclear matter is most tightly bound at about 1/10 the nuclear saturation density, somewhat below 
the density for which the pairing fields reach their maximum values. For larger and smaller values of the density, the 
pair is larger in extent and less bound. 


FINAL REMARKS 

Our numerical results indicate that the two principal physical ingredients that fix the form and magnitude of the 
pairing fields and the gap function are the energy of the nucleon-nucleon virtual state and the momentum dependence 
of its vertex function or, equivalently, the two-nucleon singlet scattering length and the momentum dependence of 
the half-on-shell zero-energy singlet T-matrix. The short-range two-nucleon correlations contained in the momentum 
dependence of the pairing fields and the virtual state are determined by the interplay between the short-range repulsion 
and medium-range attraction of the nucleon-nucleon interaction. If the short-range correlations are held fixed, the 
magnitude of the pairing fields and the position of the virtual state depend only on the overall strength of the 
interaction. 

Pairing is also expected to occur in the 3 Si- 3 Di deuteron channel. This possibility has been investigated by two 
groups Jll]{ri||3(| , which both obtained nonrelativistic S-D gap functions that reach values of about 10 MeV. Although 
we have not performed such calculations, we note that the large values obtained for the S-D gap function are not 
surprising, given the correlation we observed in Fig. [|. Since the deuteron is a bound pair while the 1 So two-nucleon 
state is virtual, the stronger effective two-nucleon potential in the 3 Si- 3 Di channel would also be expected to produce 
a much larger gap function and pairing fields. 

Several recent works have pointed out that medium polarization effects result in an important renormalization of 
the effective (non-relativistic) nucleon-nucleon interaction Q ^6| . One of the observed effects of this renormalization 
is a large reduction in the magnitude of the gap function. Such a reduction is consistent with our results if the 
polarization effects of the medium reduce the attraction of the effective two-nucleon potential. Another recent study 
has shown that the quark substructure of the nucleons and exchanged mesons is also expected to lead to a less 
attractive two-nucleon potential (37|. We plan to extend our study of pairing to next take into account such effects. 
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FIG. 1. Pairing fields and gap functions at the Fermi momentum, obtained with the parameters of Ref. [19] at two values 
of the cutoff parameter, A = 2 GeV/c (solid lines) and A = 1 GeV/c (dashed lines), for the <t-uj model. The fields are Ao, A.,, 
and A g , are in decreasing order. 


FIG. 2. Pairing fields and gap functions at the Fermi momentum, obtained with a cutoff parameter of A = 2 GeV/c using 
the parameters of Ref. [18] (dashed lines) and Ref. [19] (solid lines), for the a-co model. The fields Ao, A s , and A g , are in 
decreasing order. 


FIG. 3. Normalized pairing fields and gap functions, as a function of the momentum, for several values of the Fermi 
momentum, obtained with the a-u> model parameters of Ref. [19] and a cutoff parameter of A = 2 GeV/c. The ratio 
fio(fc,0;0)/fio(0,0;0) of the Bonn-B 1 So half-on-shell T-matrix at zero energy is also shown (short-dashed line). 


FIG. 4. Normalized pairing fields and gap functions, as a function of the momentum, at a Fermi momentum of = 0.5 
fin -1 , obtained with the parameters of Ref. [18] at two values of the cutoff parameter, A = 2 GeV/c (solid lines) and A = 1 
GeV/c (dashed lines), and the parameters of Ref. [19] with a cutoff of A = 1 GeV/c (dotted line), for the a-u> model. The ratio 
fio(fc,0;0)/tio(0,0;0) of the Bonn-B 1 So half-on-shell T-matrix at zero energy is also shown (short-dashed line). 


FIG. 5. The gap function A 9 at a Fermi momentum of Icf = 0.5 fm _1 plotted versus the momentum of the bound or virtual 
state, —ikb, for each set of interaction parameters and cutoff. Calculations with cr and u> mesons and a, u>, n and p mesons are 
labelled in the figure. The zero-range calculation is labelled ZR. 


FIG. 6. The pairing fields Ao and As at a Fermi momentum of k,F = 0.5 fm -1 plotted versus the momentum of the bound 
or virtual state, —ikb, for each set of interaction parameters and cutoff. Calculations with cr and ui mesons and cr, u>, n and p 
mesons are labelled in the figure. The zero-range calculation is labelled ZR. 


FIG. 7. The inverse of the coherence length l/£ at a Fermi momentum of kF = 0.5 fm -1 plotted versus the momentum of 
the bound or virtual state, —ikb, for each set of interaction parameters and cutoff. Calculations with cr and w mesons and a, 
u>, 7r and p mesons are labelled in the figure. The zero-range calculation is labelled ZR. 
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FIG. 8. Gap functions A g obtained for several sets of a-ui interaction parameters that have been constrained so as to yield 
a virtual state close to the physical one. The non-relativistic calculation of Ref. [30], which uses a Reid soft-core potential, is 
also shown. 

FIG. 9. Coherence lengths £ obtained for the same sets of a-uj interaction parameters constrained so as to yield a virtual 
state close to the physical one as in Fig. 8. 
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